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We show that in a turbulent flow transporting suspended sediment, the unsaturated 
sediment flux q(x,t) can be described by a first-order relaxation equation. From a mode 
analysis of the advection-diffusion equation for the particle concentration, the relaxation 
length and time scales of the dominant mode are shown to be the deposition length 
HU/Viaxi and deposition time H/V{ a \\, where H is the flow depth, U the mean flow 
velocity and Vf a ii the sediment settling velocity. This result is expected to be particularly 
relevant for the case of sediment transport in slowly varying flows, where the flux is never 
far from saturation. Predictions are shown to be in quantitative agreement with flume 
experiments, for both net erosion and net deposition situations. 



1. Introduction 

Suspension is an important mode for the transport of sediments by fluid flows. It 
occurs when the falling velocity of the particles is smaller than the turbulent velocity 
fluctuations, so that particles can remain suspended for a long time, trapped by turbulent 
eddies, before they eventually fall back on the bed due to gravity. In nature, one observes 
suspension in large rivers, i.e. in their downstream part, where large amount of fine 
particles have been collected from the catchment basin. Rivers that ordinarily present 
bed-load transport (the moving particles remain close to the bed) can also experience 
suspension (the particles are present over the whole flow depth) when the water discharge 
is unusually large, e.g. during flood events. 

Vertical concentration profiles and overall sediment fluxes are among the major issues 
- see the pioneering works of Rouse (1936) or Vanoni (1946). From the point of view 
of hydraulic engineering, the problem is satisfactorily solved for rivers in a steady state, 
although some questions are still open, such as particle trapping by turbulent eddies, 
or the structure of the flow near the bottom where the concentration is large (Nielsen 
1992; Nezu 2005). However, the response of the sediment flux to temporal or spatial 
changes of the flow is largely unknown. Such changes may be induced, for instance, 
by long gravity waves, or a sudden increase of the flow rate, or variations of the river 
slope or geometry. Two typical problems of relaxation downstream a change in the flow 
conditions are depicted in Figure [Tl which will be studied in Sj4j that of a small change 
of the slope of the bottom (Fig. Ilk), and that of a change of the bottom conditions, 
from non-erodible to erodible (Fig. fib) . The suspended sediment response is expected 



2 



P. Claudin, F. Charru and B. Andreotti 



(a) 



TTT 



flow 



slope 5 slope S + 8S 



(b) 



clear fluid 



flow 



non erodible erodible 



Figure 1. Sketches of the situations studied in 94] (a), flow over a slope change; (b), flow over 
the passage from non erodible to erodible bed. JTand S are the river depth and bottom slope 
for x < 0, and <&o(z) is the corresponding saturated concentration profile; 8H and SS are the 
variations in H and S downstream of the change at x — 0, and & sa t(z) is the new saturated 
concentration profile reached at the distance x w L sat . 



to have a strong effect on the dynamics of the erodible bottom, especially on the for- 
mation of dunes or bars, or, at larger scale, on the development of meanders (Seminara 
2006). Specific relaxation problems have been investigated by numerical integration of 
the Reynolds-averaged Navier-Stokes equations, using mixing length or k — e turbulence 
models (Hjelmfelt & Lenau 1970; Jobson & Sayre 1970; Apmann and Rumer 1970; van 
Rijn 1986a; Celik & Rodi 1988; Ouillon & Le Guennec 1996). 

In the case where bed-load is dominant, or in the aeolian situation (saltation), it has 
been shown that the evolution of the sediment flux q can be accounted for by a relaxation 
equation of the form 

T sst d t q + L sat d x q = q sst - q, (1.1) 

where q sa t is the saturated flux, T sa ± and L sat are the relaxation time and length scales 
- i.e. the time and length over which the flux relaxes toward saturation. This satura- 
tion corresponds to the homogeneous and steady state for which sediment transport is 
constant both in space and time for given flow conditions. Such a first order equation 
was first introduced in the aeolian context, as the simplest equation describing relax- 
ation effects (Saucrmann, Kroy & Herrmann 2001; Andreotti, Claudin & Douady 2002; 
Kroy, Sauermann & Herrmann 2002; Andreotti 2004). It was then shown that, for water 
flows, this equation can be derived from analysis of the erosion and deposition rates, 
the relaxation scales being there related to particle deposition (and not particle inertia) 
(Charru 2006; Lajeunesse, Malverti & Charru 2010). The importance of the relaxation 
length L sat appeared to be crucial in particular for stability analyses of the erodible bot- 
tom, for selection of the ripple wavelength (Fourriere, Claudin & Andreotti 2010). The 
importance of relaxation phenomena for suspensions in turbulent flow is well-known in 
the context of hydraulic engineering (Yalin & Finlaysen 1973; van Rijn 1986b; Celik & 
Rodi 1988; Ouillon & Le Guennec 1996). This importance has also been recognized in 
the context of geomorphology; in particular, Lague & Davy (2009) proposed a deposition 
length of sediment as the relevant transport length. However, a derivation of a relaxation 
equation of the form (1.1 1, from firm hydrodynamic grounds, is still lacking. 

In this paper, we discuss the conditions under which an equation of the form ( |1.1| 
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can be derived for turbulent flows, when suspension is the dominant mode of transport, 
with particular emphasis on the identification of the saturation length and time scales. 
The article is organised as follows. In the next Section, we present the flow models 
and saturation conditions. In Sj3]we perform a mode analysis of the advection-diffusion 
equation for the particle concentration applied to unsaturated cases, and then identify the 
saturation length and time scales of the flux. The relevance of this approach is illustrated 
in 5j4]by treating few examples: (i) the effect on the sediment flux of a change in the river 
slope; (ii) the change from a fixed to an erodible bed and (iii) the deposition of sediments 
from a source near the free surface. For the last two situations, the predictions of the 
model are tested against experimental data from the literature. 



2. Flow models 

2.1. Logarithmic flow model 

We consider the free-surface, turbulent flow of a fluid layer of thickness H over an erodible 
bed. For the sake of simplicity, we restrict the discussion to flows invariant in the spanwise 
direction, i.e. two-dimensional, with streamwise coordinate x and upwards transverse 
coordinate z. Measurements have shown that the profile of the streamwise velocity is 
close to the logarithmic law 

«.(*) = U - In ( ^ , (2.1) 

K \ Z J 

where u* is the friction velocity, zq is the hydrodynamical bed roughness, and n = 0.4 
is the von Karman coefficient. For steady flow where the shear stress is balanced by 
the streamwise component of gravity, the shear stress increases linearly from zero at the 



free surface to = pu\ at the bottom, so that the logarithmic velocity profile ( |2.1 ) 
corresponds to a parabolic eddy viscosity v t (Nezu & Rodi 1986), given by 

+ fi-JL). ( 2 .2) 



u*H V H / V H 



From (2.1 ), the depth-averaged velocity U is given by 

A s - = ± L (-) - ll , (2.3) 

U* K I \Zq J 

with typical value A = 10, corresponding to z /H 0.01 (Raudkivi 1998). 

We assume that the sediment concentration <f> is governed by the advection-diffusion 
equation 

dcj) d<t> d ( n d6\ d / dc/> ,„ \ 

it + u *£ = & ( D £ ) + oz [ D Tz + ) ' (2 - 4) 

where D is the particle eddy diffusivity and Vfaii the settling velocity. Measurements 
have shown that D(z) is reasonably parabolic and proportional to the eddy viscosity 



(2.2), with turbulent Schmidt number 



Sc=^ (2.5) 

in the range 0.5-1 (Coleman 1970; Celik & Rodi 1988; Nielsen 1992). The settling 
velocity Vf a n is taken uniform, and, when needed for comparison with experiments, equal 
to that of a single particle in quiescent fluid. Note that the above modelling ignores 
inertial effects on particle motion, in particular their ejection from the core of vortices 
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and their clustering (Bee et al. 2007; Hunt et al. 2007). We also limit the discussion 
to dilute suspensions, i.e. small volumic particle concentration (f>, for which there is no 
significant feedback of the particles on transport. 



Solving (2.4 1 requires two boundary conditions, one at the free surface and one on the 



sedimentary bed. At the free surface, the net vertical flux vanishes, giving 

M 



D 



dz 



4>v f 



fall 







at 



H. 



(2.6) 



At the bottom, just above the bedload layer where particles mainly roll and slide on 
each other, the diffusive flux is equal to the erosion flux tp^ } i.e. the volume of particles 
entrained in suspension per unit time and bed area (Parker 1978; van Rijn 1986a): 

M 



D 7T = 
az 



fl- 



at z = 0. 



(2.7) 



The erosion rate, or 'pickup function', is generically an increasing function of the basal 
shear stress above a threshold. Its functional form is determined phenomenologically 
from experiments and depends on the nature of the bed - e.g. whether it is consoli- 
dated/cohesive or not, composed of grains or containing clay, etc. (Shields 1936; Einstein 
1950; Engelund 1970; van Rijn 1984b; Hanson & Simon 2001; Briaud & al. 2001; Bonelli 
et al. 2007). 

Two remarks have to be made here. First, the bottom condition (2.7) applies for 



steady and homogeneous as well as unsteady or heterogeneous flows. In the latter case, 
the erosion flux may be different from the deposition flux, so that the net flux is nonzero, 
which may lead to variations of the bed topography (not necessarily, as in the experiments 
to be discussed later). Possible variations in the bed topography will be ignored here. 



Second remark, the boundary condition (2.7 1 corresponds to a bed allowing unlimited 



sediment supply. For more general situations (e.g. fixed bed), slightly different boundary 
conditions have been proposed, see Celik & Rodi (1988), which however requires an 
empirical constant or reference concentration near the bottom to be given, which varies 
along the channel. Finally, assuming that particles have the same mean velocity as the 
fluid, u x , the flux of suspended particles, per unit length in the spanwise direction, is 
given by 

q= { (pu x dz. (2.8) 



The concentration equation ( 2.4 ) with the boundary condition ( 2.6 1 admit a steady and 



homogeneous solution corresponding to the balance of the settling and diffusive fluxes, 

z/H 



$sat(z) = $b 



1 



1 + Z/ZQ 



1 + zo/H 



where /3, known as the Rouse number, is defined as /3 
concentration $b = ^sat(O) is determined from the condition (2.7) as 



(2.9) 

(ScVfaii)/(«M*) and the bottom 



$h = 



all 



(2.10) 



Note that (2.9) differs slightly for the classical expression of the Rouse profile (Nielsen 
1992) because the location where the velocity ( |2.1| vanishes and where the bottom bound- 
ary condition (2.7) applies has been chosen to be z = instead of z — zq. Suspension 
typically occurs when Vfaii < 0.8 u* (Freds0e & Daigaard 1992), which corresponds to 
/3 < 4/3 with Sc = 2/3. Figure [2^i displays the velocity profile (2.1), normalized by 
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Figure 2. (a) Logarit hmic flow model: normalized velocity profile (dashed line) and normalized 
concentration profiles (|2.9| for three values of /3 (solid lines); (b) plug flow model: normalized 
concentration profiles Q2.12[| for the corresponding values of a = 6/3. 



u x (H), and the concentration profile (2.9), normalized by $b, for three typical values of 
P. 

2.2. A simplified plug flow model 

In order to get analytical results, a simplified plug flow model will be used in the following, 
which appears to provide accurate results as long as assessment of the relaxation equation 
(1.1 ) is pursued. This model corresponds to uniform flow velocity u x {z) — U and friction 
velocity = U/X, where A is the same constant as in the previous section. Such a plug 
flow model is of course a rough description, and the velocity profile actually does not 
switch from zero on the bed to its average value on a vanishing vertical distance, leading 
to an infinite shear. This problem is not present in the logarithmic model, which is more 
realistic from this point of view. However, as shown below, these two models do not differ 
much as far as the relaxation modes are concerned, which means that what occurs very 
close to the bed is not very important for the present purpose. Accordingly, a uniform 



particle diffusivity Dq will be taken in the concentration equation (2.4 1 and boundary 
conditions ( 2.6|2.7 ), equal to the average of the parabolic distribution given by (2.2 1 and 
(2.5). Up to a small correction of order zq/H, the diffusivity Dq is given by 

Do n 
u*H 6Sc 



(2.11) 



With this uniform diffusivity Dq, the advection-diffusion equation (2.4) admits the 
steady and homogeneous solution 



3>b exp 



(-4) 



with a = 



a = — = 6/3, 



(2.12) 



which also satisfies the boundary condition (2.6) at the free surface. The boundary 
condition at the bed (2.7) determines the bed concentration (2.10). Figure [2|d displays 
the concentration profile (2.12), normalized by $b, for three typical values of a = 6/3. It 



can be seen that for the same value of the Rouse number, the plug flow model predicts 
sediment concentration slightly larger than that of the logarithmic flow model. Small 
values of a correspond to strong suspensions, i.e. situations for which the sediment is 
distributed almost uniformly over the whole depth of the flow. This is achieved when 
the settling velocity is small (very fine particles) or when the diffusivity is large (large 
flow velocity). Finally, the saturated particle flux per unit width q sat , normalized by the 
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water flux UH, is given by 



For small a, this dimensionless flux tends to as expected. 



3. Non- homogeneous and unsteady flows 



In this section, we successively consider a spatial evolution problem (1 3.1 1 and a tern 



poral evolution problem ((3.2). These problems are solved using a mode analysis, i.e. 
the departure of the concentation field from the saturated distribution $> sa t(z) is de- 
composed as a sum of terms of the form c(x,t)f(z). It is shown that, for the spatial 
problem, there exists a discrete set of amplitudes c n {x) oc e~ x ^ Ln , and for the temporal 
problem, there exists a similar set of amplitudes c n {t) cx e _t / T ™. Then the sediment flux 
is shown to be dominated by the mode with the largest length or time, L\ or T\. This 



result demonstrates that for large scale problems, the relaxation equation (1.1) retains 
the most important features of unsaturated sediment transport, with relaxation scales 
L sa t and T sat equal to the largest scales arising from the mode analysis. 

The analytical calculations presented below use the plug flow model because of its 
simplicity, in the spirit of the work of Mei (1979). Calculations for the logarithmic model 
are not reported in detail, but the corresponding results are plotted for comparison in 
some of the figures. 

3.1. Spatial evolution and the relaxation lengths 

We consider the situation where, for given flow conditions and corresponding saturated 
concentration profile $ S at(z), the actual concentration profile at some point, say x = 
0, is <&sat(z) — 4>{ x — 0, z) where <j>(x,z) is a 'concentration defect'. We search for 
the distance at which the saturated distribution < I > sa t(z) is recovered, corresponding to 
vanishing 4>{x, z). Looking for normal modes of relaxation of the concentration defect of 
the form 

^,z)=$ b /(z)exp(-.T/L), (3.1) 



we get from the advection-diffusion equation (2.4) 



Au» D\ , d ( df 



L < Li) f + T Z { D t +fV ™)=°- (3 ' 2) 



At the free surface, the zero flux condition (2.6) gives 



£>~+/VWl = at z = H. (3.3) 

dz 

On the bed the friction velocity u* is assumed to be uniform. The erosion flux ip^, which 
depends only on , is uniform too. Hence, the disturbance of (p+ is zero, so that, from 



= at z = 0. (3.4) 

dz 

The above differential problem is solved numerically for parabolic D (logarithmic flow 
model) and analytically for uniform D — Dq (plug flow model). For uniform Dq, equation 



(3.2) has solutions of the form f(z) oc exp(Kz/H), where K has to satisfy a quadratic 
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Figure 3. (a) Variation with a of the three smallest roots of (3.6 1, (b) Corresponding 
relaxation normalised lengths L n /Ld for X/K, = 50; plug flow ( — ); logarithmic flow ( 



equation with roots K + and K_ given by 



01 . , t , X H H 2 a 2 
K+= ±iKi with Ki = \ 

2 1 1 V K L L 2 4 



(3.5) 



Then the boundary conditions ( 3.3||3.4 ) select a discrete set of relaxation lengths L ver- 
ifying : 

tani^d-^)" 1 . (3.6) 

This equation has an infinite number of real positive solutions K- m , n > 1. Figure [3^, 
shows the variation with a of the three smallest ones (n = 1,2,3). For small a, these 
solutions behave as K\\ ~ ^/a and ifi„ ~ (n — l)ir for n > 2. 



The corresponding relaxation lengths L n are found from (3.5) 

A 



H _ 1 

L 7r 2 




■4K?. 



n > 1. 



(3.7) 



They are displayed for n = 1,2,3 in Figure [3p as a function of a (solid lines), normalised 
with the characteristic deposition length 



V'faU 



(3.8) 



It can be seen that L\ is much larger than the higher-order relaxation lengths - typically 
by one order of magnitude. Remarkably, in the limit of small a (large flow velocity or 
small settling velocity), the largest length L\ tends to Ld, whereas higher-order lengths 
remain on the order of the flow depth H: 



X/K 



(n 



: H for n > 2 



(3.9) 



Figure ^jp also displays the normalized relaxation lengths obtained from the logarithmic 
flow model (dashed lines), from numerical integration of (3.2 1 with parabolic D. It can 



be seen that these lengths are close to those from the plug flow model, especially for the 
largest length L%. Note that L n /Ld is weakly sensitive to the value of X/IC: doubling 
this ratio does not bring any visible change, at least for a < 2. 
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The eigenfunctions f n {z) are given by 



fn(z) 



cos ( K\ 



H 



2K„ 



$\n\K in ^ exp( 



a z 
2 H 



n > 1, 



(3.10) 



with the normalization condition /„(0) = 1. These eigenfunctions are displayed in Fig- 
ure[4]for n = 1, 2, 3 (solid lines), for a = 0.1 (Fig. [4^,) and a = 1 (Fig.|4j3). It can be seen 
that fi(z) decreases slightly and monotically from bottom to top, whereas higher-order 
eigenfunctions oscillate, more and more strongly with increasing n. Figure [4] also displays 
the eigenfunctions from the logarithmic flow model (dashed lines). It can be seen that for 
the mode associated with the largest length L\ (n = 1), eigenfunctions of both models 
remain very close to each other, and that differences become larger as n increases. 

Let us turn to the sediment flux. The contribution of the n th -eigenmode to the sed- 
iment flux, Q„, normalized with the characteristic sediment flux UH&b and the expo- 
nential x-dependence, is 

cH 

f n {z)dz. (3.11) 



1 



1 _ 1 

$ h UH ® n ~ H 



exp(—x/L r , 

Table [l^, displays the contribution of each of the first three modes to the sediment flux, 
i.e. the right-hand side of the above equation. It can be seen that the contribution 
of the first mode n = 1 strongly dominates. The smallness of the contribution of the 
higher-order modes is due to the oscillations of the eigenfunctions, as shown in Figure |4j 
For small a, the normalised flux is close to one for n = 1 and decreases as a/((n — l)ir) 2 
for n > 2. 

The general form of the concentration defect finally is 



x,z) = $ b V a n f n (z)exp(-x/L n ), 



(3.12) 



where the relaxation lengths L n are given by ( 3.6|3.7 l, the eigenfunctions f n (z) are given 
by (3.10), and the coefficients a n have to be determined by the concentration profile 
imposed at x = 0. Such a determination will be illustrated in section |4j 

3.2. Temporal evolution and relaxation times 

We now consider an unsaturated concentration profile at initial time t — 0, say & S a,t(z) — 
(f>(z,t = 0), uniform in the strcamwise ir-direction, and search for the time needed for 
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n 


1 


2 


3 






A 1 


M 


A 3 


^4 


(a) 


a = 0.1 
a = 1 


0.9836 
0.8533 


0.0099 
0.0832 


0.0025 
0.0240 


(b) 


a = 0.1 
a = 1 


0.967 
0.724 


0.020 
0.150 


0.005 
0.047 


0.002 
0.022 



Table 1. (a) Contrib ution of the three lowest-order eigenmodes to the normalized sediment 
flux (r.h.s. of equ ation (3.11 1), for a = 0.1 and a = 1. (b) Normalized coefficients A n = a„/5$b 
of the expansion (4.2|, computed from a projection over four modes. 



relaxation to the saturated distribution $ sa t(^) given by (2.12), i.e. vanishing concentra- 

ection, 

(3.13) 



tion defect <p(z,t). Calculations go along the same lines as in the previous sub-section, 
so they are only briefly sketched here. Looking for normal modes of the form 



<Hi, z) = i b g{z)exp(-t/T), 



we get from the equation (2.4) the equation governing the eigenfunctions g{z): 



T 



9 - 



d 

dz 



A) ~r + g vfaii 



0. 



(3.14) 



This equation has solutions of the form g(z) oc exp(Kz/H), where K has to satisfy a 
quadratic equation with roots K + and if_ defined as 



a 



± LfTi with ifi 



JCu^T 



T 



(3.15) 



The boundary conditions at z = and z = H are the same as in the previous section, so 
that A'; still verifies equation (3.6), with same solutions K mi n > 1. The corresponding 

by 

2 

= i?f?+^, u>l. (3.16) 



relaxation times T n are then given by 

H 



Introducing the characteristic deposition time 

d -F fall - U ' 

the relaxation times are, in the limit of small a, 

aid in 



Ti ~ T d , 



for n > 2. 



(3.17) 



(3.18) 



(n- 1) 2 tt 2 

As for the spatial problem, the sediment dynamics is dominated by the largest time T\ , 
equal to the deposition time Td for strong suspensions. 



4. Two illustrations, and comparison to experiments 

4.1. Effect of a change in the bed slope 

Consider the situation depicted in Figure [T^,, of a flow with saturated concentration 
profile &o(z) which experiences a small variation SS in the bottom slope at x = 0, either 
positive or negative. This variation leads to a small change of the water depth and 
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0.5 1 0.5 1 



Figure 5. Relaxation of the concentration defect <j>(x,z) after a slope change as sketched in 
Figure nk, with summation over the three first modes (plug flow model), (a), Profiles of <p(x, z) 
at the four downstream positions x/Ld = 0, 0.3, 1 and 3, for a — 0.1; (b), same for a = 1. 



friction velocity, according to 8H/H = —5u*/u* = —^SS/S. This change occurs on a 
hydrodynamic lengthscale Lh given by the balance between the acceleration USU/Lh 
and the force gSS, i.e. L^/Ld = UVf a \i/2u1 = XK.a/2. The present analysis is valid for 
small Lh/Ld, a condition which is fulfilled for small a. 

The change in the saturated concentration profile due to the slope variation SS is, at 
the linear order in 5u*. 



5$(z) = cxp (-a-^j with 5®k = 



VfaXl 

where </5^(u«) is the derivative of the erosion rate ip^(u*). The concentration defect at 
x = corresponds to this change (<j)(0,z) — 6<&(z)), so that the coefficients a„ of the 



expansion (3.12| must satisfy 



S<i> 



OO 

b cxp (-^Jj) =Y1 a ™ cos ( 



n=l 



Z 



2K U 



sm 



(4.2) 



These coefficients can be determined from the projection of the above equation on the 
eigenfunctions, i.e. truncation of the sum on the r.h.s. to p terms, multiplication by 
each eigenfunction, and integration of both sides from to H. A linear system of p 
equations is obtained, whose solution gives the coefficients ax, — , a p . Table [T]} displays 
the normalized coefficients A n = a„/(5$b resulting from the projection over p = 4 modes, 
for two values of a. One can see that the first mode captures most of the weight; the 
contribution of the second one is smaller, but still significant, and higher modes are 
negligible. We have checked that considering more terms in the expansion has negligible 
effect on the dominant coefficients. The small weight of the oscillating modes is consistent 



with the slow variation with z of the initial concentration profile (4.1). 

Figure [5] displays profiles of the concentration defect <p(x, z) at the location of the 
slope change, x/Ld = 0, and three positions downtream, for a = 0.1 (Fig. [5^l) and 
a = 1 (Fig.|5p). It can be seen that at the position x/Ld — 3, the concentration defect is 
nearly zero. The flux defect, i.e. the depth-integrated profiles of the concentration defect, 
correspondingly decays towards zero (not shown) and does so almost exponentially with 



relaxation length Ld, as predicted by equation (1.1). This confirms that the dominant 
mode with relaxation length Ld captures most of the sediment flux variations. For this 
example, as well as for the next ones, the question of the plug flow limit is not crucial: 
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| van Rijn (1986b) | Ashida & Okabe (1982) | Jobson & Sayre (1970) | 







run 5 






run 6 


runs FS1 and FS1A 




erosion 


erosion 




deposition 


deposition 


Symbol 


(*) 


(°) 






(□) 


(■) 


H (cm) 


25 




4.3 




40.7 


U (cm/s) 


67 




37.3 




29.1 


tt, (cm/s) 


4.77 




3.63 




4.48 


Vfaii (cm/s) 


2.2 




1.85 




1.0 - 2.0 


L d (m) 


7.6 




0.87 




5.9 - 11.9 


a 




3.1 






2.5 


1.2 


Sc 




0.41 






0.33 


0.18 - 0.36 



Table 2. Hydraulic parameters H, U, u, and Vf a ii of the experiments, and Ld = ((7/Vf a n)//, a 
from the exponential fit of the downstream concentration profile, and Sc = /ta«*/6Vf a ii. 



as far as the relaxation modes are concerned, the logarithmic and plug models do not 
differ much (Fig.|4|. 



4.2. Net erosion experiments 

Another situation of interest is that of a flow of clear fluid on a non-erodible bed (§o(z) = 
0) reaching an erodible bed lying in x > 0, as sketched in Figure [T|d. Suspension develops 
downstream until the saturated concentration profile $ sa t( z ) is reached. The analysis 
goes along the same lines as for the slope change. The concentration defect <fi(x,z), 
can be decomposed on the eigenfunctions (3.10), with coefficients a n determined by the 
concentration profile at x = 0. The equation to be satisfied turns out to be the same as 
(4.2) with <£>b instead of <5$b- Thus the coefficients are a n = A„(S$b with the normalised 
coefficients A n given in Table [Tja. 

The prediction that the eigenmode with the largest relaxation length captures most 
of the sediment flux can be assessed from the experimental observations of van Rijn 
(1986b) and Ashida & Okabe (1982) - non-Japanese readers can access these latter data 
in the paper of Celik & Rodi (1988). These experiments precisely correspond to the 
sketch depicted in Figure [TJb. Their hydraulic parameters are given in Table [2j The 
spatial evolution of the concentration profiles has been measured at different locations 
downstream the transition point at x = 0. The corresponding sediment flux q, which 
is zero for x <C 0, increases downstream until it reaches the saturated value Q S at- 

We 

determined this flux from integration of the measured concentration profile at each al- 
location. From the hydraulic parameters, the deposition length can be computed as 
Ld = {U /Vfaw)!! - in the following we will not distinguish between L\ and Ld, although 
they can differ by sa 20% for a on the order of unity (Fig.^). It appeared that for Ashida 
& Okabe (1982) the location of the farthest downstream measurements corresponds to 
x/Ld — 8.1, which is large enough for the sediment flux to be saturated. Figure [6^, 
displays the corresponding concentration profile, and an exponential fit providing, from 
(2.12 1, the value of the parameter a reported in Table[2] For van Rijn (1986b), we found 
x/Ld = 1-3, not large, preventing any straightforward determination of the parameter a. 
For both experiments, the saturated flux q sat was estimated as that providing the best 
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Figure 6. Net erosion experiments. (a) Concentration profile at the farthest location 
x/Ld = 8.1 measured by Ashida & Okabe ( 1982, run 5), normalized with the depth-averaged 
concentration $ re f; solid line: $sat given by ( |2.12[) with a — 3.2. (b) Relaxation to saturation 
of the normalized sediment flux versus x/Ld; symbols: van Rijn (1986b) and Ashida & Okabe 
(1982) (see Table 12} ; solid line: exponential relaxation ( |4.3[ ). 

fit to the exponential curve 

^- = 1 - exp(-x/L d ). (4.3) 

<7sat 

Figure^ diplays the variation of q/q sa .t with x/Ld] it can be seen that the data points 
fall quite well on the exponential curve. Note that the deposition lengths Ld differ by one 
order of magnitude between the two experiments. The data collapse therefore supports 
a first-order relaxation process with characteristic length equal to Ld- 

4.3. Net deposition experiments 

Ashida & Okabe (1982) have also performed experiments in which the initial concentra- 
tion profile is oversaturated (run 6), i.e. the initial sediment flux q at x = is larger 
than g sa t , so that the sediment settle until the saturated regime is reached further down- 
stream. Figure [TJd displays the sediment flux, obtained from the measured concentration 
profiles, together with the exponential relaxation curve now given by 

— = 1 + ( — ~ A cxp(-x/L d ) (4.4) 

<7sat \ <Zsat / 

where qo and g sa t were determined by curve fitting, and Ld is given in Table [2j Again, 
the agreement is quite good, showing that the mode with relaxation length Ld captures 
most of the deposition process. 

As a confirmation, Figure [7^, compares the concentration profile measured at the lo- 
cation x = to its projection over one single mode. This projection is computed from 



an empirical representation of this inital profile, using the expansion (3.12) and the sat 



urated flux <& sat (z) measured from the concentration profile at x = 8.1 Ld, fitted by the 



exponential form (2.12) with a = 2.5. It can be seen that the resulting profile is in good 
agreement with the measurements. Note that a — 2.5 corresponds to Schmidt number 
Sc = anu*/6Vi a \\ — 0.33, which is slightly below the usual range 0.5-1 (Coleman 1970; 
Celik & Rodi 1988; Nielsen 1992). 

Other net deposition experiments have been performed by Jobson & Sayre (1970). In 
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Figure 7. Net deposition experiments of Ashida & Okabe (1982) (run 6, see Table[5f . (a) Con- 
centration profiles at x = normalized with the depth-averaged concentration $ rc f; solid line: 
$o reconstructed from $ sat and one single mode for the concentration defect, (b) Normalized 



sediment flux versus x/Ld, experiments and exponential relaxation (4.4 1 



this work, the particles were released near the water surface, so that the initial concen- 
tration profiles exhibits a peak close to z = H, as shown in Figure [8^i. In contrast to 
Ashida & Okabe experiments, expanding the concentration defect over one single mode 
is not sufficient to get a good representation of this profile; an expansion over four modes 
provide a much better description, as shown in Figure [8^,. However, the high-order modes 
are expected to vanish over a short distance, on the order of a few flow depths H, and 
the exponential relaxation to be recovered at large distances. This scenario is evidenced 
in Figure [8Jd, which displays the normalized flux versus x/H, measurements and the 
exponential curve (4.4). Here, due to uncertainties on the falling velocity (see Table [2] of 
the present paper and Figure 6a of Jobson & Sayre 1970), the deposition length L d was 
determined, together with qo and g sa t, by fitting the experimental data points. We found 
Ld = 5.2 m, which is close to the range of the expected values displayed in Table [2j 
although slightly smaller. We finally note that in the course of the reconstruction of 
$0(2), we found a = 1.2 from the saturated concentration profile, which corresponds to 
Schmidt number in the range 0.18-0.36, slightly smaller, again, than the usual range. 



5. Concluding remarks 

In this paper, we have discussed the conditions under which a first-order relaxation 
equation for the sediment flux q can be derived for turbulent flows, when suspension is 
the dominant mode of transport. From a mode analysis of the linear advection-diffusion 
equation for the particle concentration, it was shown that the sediment flux is dominated 
by the mode corresponding to the largest relaxation length for spatially varying flows, 
or the largest relaxation time for time-dependent flows. These relaxation scales were 
identified as the deposition length i?f//Vf a ii and the deposition time H/Vf a \\, where H is 
the flow depth, U the mean flow velocity and Vf a ii the sediment settling velocity. This 
result is expected to be particularly relevant for the case of sediment transport in slowly 
varying flows, for which the flux is never far from saturation. Predictions of the sediment 
flux were shown to be in quantitative agreement with flume experiments, for both net 
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Figure 8. Net deposition experiments of Jobson & Sayre (1970) (see Table EL (a) Measured 
concentration profile at x — normalized with the depth-averaged concentration $ rc f; dashed 
and solid lines: <J>o reconstructed with one single mode and four modes, respectively, and a = 1.2. 
(b) Relaxat ion to saturation of the normalized sediment flux versus x/Ld\ solid line: exponential 
relaxation ( |4.4| ) with = 5.2 m. 

erosion and net deposition situations, and deposition lengths spanning over one order of 
magnitude. 

As discussed in the introduction, the relaxation equation ( |1.1| allows for the description 
of both bed load and suspended load. However, these modes of transport correspond to 
very different physical lengthscales. For bed load, the relaxation length L sat is on the 
order of 10 grain diameters (Fourriere, Claudin & Andreotti 2010). As soon as the 
flow depth H is larger than a few L sat) the - unstable - flat bed is insensitive to the 
presence of the free surface, and current ripples emerge at a centimetric wavelength 
(« 10L sa t). When suspended load is the dominant type of transport, we have shown 
that the relaxation length L sat is on the order of 10-100 H, which is typically four to 
five orders of magnitude larger than for bed load. Suspended transport thus prevents 
the formation of bedforms with wavelength smaller than H , and patterns such as bars, 
antidunes and meanders can be expected to emerge from linear instability, with large 
wavelengths on the order of 100-1000 H . Further work is required for the experimental 
and theoretical investigations of these instabilities. 
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